Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 14 January 2013 



(MN M£X style file v2.2) 



Asymptotic Giant Branch Stars in the Sculptor Dwarf 
Spheroidal Galaxy 



o 

(N 



John W. Menzies 1 , Michael W. Feast 2 ' 1 , Patricia A. Whitelock 1 ' 2 and Noriyuki 
Matsunaga 3 

1 South African Astronomical Observatory, P.O.Box 9, 7935 Observatory, South Africa. 

2 Astronomy, Cosmology and Gravitation Centre, Astronomy Department, University of Cape Town, 7701 Rondebosch, South Africa. 

3 Kiso Observatory, Institute of Astronomy, The University of Tokyo, 10762-30, Mitake, Kiso, Nagano 397-0101, Japan. 



OS 

O 

u 

6 



> 

m 



O 



14 January 2013 



X 



ABSTRACT 

JHKs photometry is presented for a 35 arcmin square field centred on the Sculptor 
dwarf spheroidal galaxy. With the aid of published kinematic data definite galaxy 
members are identified and the width in J — K of the colour-magnitude diagram is 
shown to be consistent with an old population of stars with a large range in metal 
abundance. We identify two Asymptotic Giant Branch variables, both carbon Miras, 
with periods of 189 and 554 days, respectively, and discuss their ages, metallicities 
and mass loss as well as their positions in the Mira period-luminosity diagram. 
There is evidence for a general period-age relation for Local Group Miras. The mass- 
loss rate for the 554-day variable, MAG29, appears to be consistent with that found 
for Miras of comparable period in other Local Group galaxies. 

Key words: stars: AGB and post-AGB - stars: carbon - stars: variables: other - 
galaxies: dwarf - galaxies: individual: Sculptor-Local Group - galaxies: stellar content. 



1 INTRODUCTION 

This paper is one of a series aimed at finding and char- 
acterizing luminous Asymptotic Giant Branch (AGB) vari- 
ables within Local Group galaxies; it follows similar work 
on Phoenix, Fornax and Leo I (Menzies et al. 2002, 2008; 
Whitelock et al. 2009; Menzies et al. 2010). 

Sculptor was the first dwarf spheroidal galaxy to be 
discovered (Shapley 1938). Variable stars were found in it 
by Baade & Hubble (1939) and Thackeray (1950) confirmed 
that these were mainly RR Lyrae variables. Establishing the 
nature of these variables was important for Baade's general 
concept of populations I and II (see letter from Baade to 
Thackeray 1950 Mar 27 in Feast 2000). Further work on the 
variables was carried out by van Agt (1978) and Kaluzny 
et al. (1995). Most of the variables are RR Lyrae stars but 
van Agt also discovered some red variables one of which we 
show below to be a Mira variable. 

Carbon stars were discovered in the system by Az- 
zopardi, Lequeux & Westerlund (1985, 1986), while Frogel 
et al. (1982) reported another three possible examples. In 
more recent times the system has been extensively studied 
for stellar chemical abundances (e.g. Coleman et al. 2005, 
Battaglia et al. 2008, Kirby et al. 2009) kinematics (e.g. 
Westfall et al. 2006; Walker et al. 2009) and star formation 
rate (Hurley-Keller 2000). 



Sculptor is believed to contain a predominantly old pop- 
ulation but with a small population of age ~ 2 Gyr and less 
(Revaz et al. 2009, especially their fig. 13). Work by vari- 
ous authors has shown that there is a large range in [Fe/H]. 
Kirby et al. (2009) found a mean value of —1.58 and an 
asymmetric tail extending out to ~ —3.0. There is evidence 
of a metallicity gradient in the galaxy with the more metal 
strong stars preferentially concentrated to the centre (see 
Tolstoy et al. 2004; Westfall et al. 2006; Walker et al. 2009; 
Kirby et al. 2009). 



2 OBSERVATIONS 

Observations were made with the SIRIUS camera on the 
Japanese-South African Infrared Survey Facility (IRSF) at 
SAAO, Sutherland. The camera produces simultaneous J, 
H and Ks images covering an approximately 7.2 x 7.2 ar- 
cmin square field (after dithering) with a scale of 0.45 arc- 
sec/pixel. Images were obtained on a 5 x 5 grid centred on the 
nominal centre of the galaxy, 01:00:06, -33:44:13 (2000.0) 
(van Agt 1978). With overlaps between adjacent fields, the 
area observed was approximately 35 arcmin squared. An ad- 
ditional field was included later to allow one of the carbon 
stars (ALW1) to be observed. 

Images were obtained at 25 epochs spread over 4 years 
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in the central 9 fields, and between 16 and 20 epochs in 
the remaining fields. In each field, 10 dithered frames were 
combined after flat-fielding and dark and sky subtraction. 
Exposure times were either 30 or 20 s for each frame, de- 
pending on the seeing and on the brightness of the sky in 
the Ks band. Photometry was performed using DoPHOT 
(Schechter et al. 1993) in 'fixed-position' mode, using the 
best-seeing if-band image in each field as templates. Aladin 
was used to correct the coordinate system on each template 
and RA and Dec were determined for each measured star. 
This allowed a cross-correlation to be made with the 2MASS 
catalogue (Cutri et al. 2003), and photometric zero points 
were determined by comparison of our photometry with that 
of 2MASS. In each field, stars in common with the 2MASS 
catalogue with photometric quality A in each colour were 
identified and the IRSF zero point was adjusted to match 
that of 2MASS. The mean number of common stars per field 
was 12 and the mean standard deviation over all fields of the 
differences between IRSF and 2MASS are 0.04 mag in J and 
H and 0.06 mag in Ks- No account was taken of possible 
colour transformations, such as in Kato et al. (2007). Those 
transformations were derived using highly reddened objects 
to define the red end and it is not obvious that the same 
transformations will apply to carbon stars. 



3 COLOUR-MAGNITUDE AND 
COLOUR-COLOUR DIAGRAMS 

Fig. [T] shows the Ks — (J — Ks) diagram and Fig. [2] the 
(J — H) — (H — Ks) diagram for all the stars discussed in 
this paper. Mean magnitudes from all epochs are used in the 
plots. 

Following the interpretation of the similar diagram for 
the Fornax dwarf spheroidal, the main features of Fig.[T]are: 
a well defined RGB with a tip near K = 13.5 mag; bluer 
stars of all magnitudes which are likely foreground stars; 
a few stars above the TRGB which are likely foreground 
stars, though some might be AGB members; two bright red 
variables; unresolved background galaxies with J — K ~ 1.5 
mag extending from K ~ 16.0 mag to fainter magnitudes. 

Bellazzini (2008) gives a relation between Mk and 
(J — K) for the tip of the red giant branch (TRGB) which 
is shown as a line in Fig. [1] assuming a distance modulus of 
19.46 mag and adjusted for consistency with an LMC dis- 
tance modulus of 18.39 mag . It is interesting to note that, 
apart from the variables, all the C stars lie below the TRGB. 
Three stars with [Fe/H]< —3.4 are shown in the figure, and 
their details are listed in Table [T] for convenience. 

We have found only two periodic red variables in Sculp- 
tor; one is our catalogue number 12002, found to be vari- 
able by van Agt (1978) and catalogued as V544. The other 
is newly discovered (our 3001), and was found to be a C 
star, designated MAG29, in a survey of 2MASS red stars by 
Mauron et al. (2004). The light curves are shown in Fig. [3] 
Following Whitelock et al. (2006) we classify them both as 
Mira variables, because they are periodic and have ampli- 
tudes AK > 0.4 mag; MAG29 has a period of 554 days, 
and V544 a period of 189 days. There is some evidence for 
long term variations in the MAG29 lightcurve, but with only 
about two cycles covered by our observations it is difficult 
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Figure 1. Colour- magnitude diagram for all stars measured in 
our Sculptor field; small grey dots are field and red giant branch 
stars; squares are variables. Black points are carbon stars. Two 
carbon stars not observed by us are plotted at their 2MASS posi- 
tions. Three stars with [Fe/H]< —3.4 are shown as blue triangles. 
The dashed line shows the approximate expected position of the 
tip of the red giant branch assuming a distance modulus of 19.46 
mag (see text). 





Figure 2. Two-colour diagram for Sculptor; symbols as in Fig.[T] 

The solid line represents the locus of Galactic carbon Miras (see 
Whitelock et al. 2009). 
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to comment constructively on this point; more data are re- 
quired. 
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Table 1. Red variables, Carbon stars and very metal-poor stars in Sculptor. 
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0.89 


ALW2 


Y 


C 


3 


01:00:10.7 


-33:40:58 


1023 


14.95 


14.33 


14.05 


0.90 


ALW4, FBMC12 


Y 


W 


3,4 


01:00:20.7 


-33:45:18 


1003 


14.60 


13.93 


13.66 


0.94 


ALW5, FBMC4 


Y 


S 


3,4 


01:00:27.4 


-33:53:04 


8008 


15.03 


14.36 


14.17 


0.87 


ALW6 






3 


00:58:27.1 


-34:01:19 




14.50 


13.77 


13.57 
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S 1020549 
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0.56 
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01:00:01.1 
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20034 
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16.30 


0.56 
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6 



Notes: 

(1) MAG29 is 2MASS00595367-3338308. It is a C star according to Mauron ct al. (2004) 
MAG29 has a peak-to-peak K amplitude of 0.87 mag. 

(2) V544 from van Agt (1978). Peak-to-peak K amplitude 0.42 mag. 

(3) ALW numbers from Azzopardi ct al. (1986). Note ALW2 and 3 are reversed in Azzopardi et al. (1985). 

(4) FBMC numbers from Frogcl ct al. (1982). 

(5) 2MASS data for ALW7 (J00582732-3401186) and ALW 8 (J00583082-3328401) outside our field. 

(6) [Fe/H]: -3.81 (Scl2030) Frebel et al. (2010); -3.48 (Scl20034), -3.96 (Scl20028) (Tafelmeyer ct al. 2010). 
Member (Y) according to: Coleman ct al. (2005): C, Schweitzer et al. (1995): S, Walker ct al. (2009): Wa, 
Westfall ct al. (2006): We. 



Table 2. JHK$ for stars in Sculptor. 



a & Scl# J H K s J-K s Memb 

(equinox 2000.0) - (mag) - 
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-33:40 
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26013 


15.15 
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C We 
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17.35 


16 


83 
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0.62 
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:58:40. 
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:44 
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17.17 


16. 


.58 


16.54 


0.63 
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16 


24 


16.14 


0.65 
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:58:43 
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-33:35 


:59 
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16 


.11 


16.03 


0.63 


Wa C 


00 
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:45 
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16 
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16.43 


0.66 


Wa 
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,7 


-33:41 


:56 
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14.66 


14 
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We 
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:49 


14018 


16.31 


15 


,74 
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0.70 


Wa We 


00' 


:58:45 


.5 


-33:38 


:02 


15026 


16.73 


16 


19 


16.16 


0.57 


Wa 


00 


:58:45 


.6 


-33:35 


:10 


15059 


17.61 


17 


.16 


17.09 


0.52 
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Notes: 

Membership references: Battaglia et al. (2008): B, Coleman et al. (2005): C, 
Kirby et al. (2009): K, Schweitzer ct al. (1995): S, Walker et al. (2009): Wa, 
Westfall ct al. (2006): We. 

This is a partial table to indicate the format and contents. The full table 
including definite members and non-members, as well as stars with unknown status 
is available as an attachment in the online version. 
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phase 

Figure 3. Kg light curves of (bottom panel) MAG29 and (top 
panel) V544 phased with periods 554 days and 189 days, respec- 
tively, and arbitrary epoch of JD 2450000. For MAG29 filled cir- 
cles (black) show first pulsation cycle, triangles (magenta) show 
second cycle, squares (cyan) show third cycle. For V544, these 
same symbols show the first, second, third pairs of cycles, respec- 
tively. 
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Figure 4. All stars observed in the Sculptor field (grey dots). 
Variables are filled squares (the northern one is V544, the one 
near the centre of the field is MAG29). Carbon stars are large 
black symbols except for two carbon stars that fall outside our 
observed field that are indicated by red filled circles on the right 
side of the figure. Three stars with very low metallicities are shown 
as blue triangles. 



Of the eight carbon stars identified by Azzopardi et al. 
(1986), six fall in the area of our survey. In a study of late- 
type stars in Sculptor, Frogel et al. (1982) find two definite 
C stars (in common with Azzopardi et al.) and three pos- 
sible ones, which however were not confirmed by Richer & 
Westerlund (1983). Our data for the variables and C stars 
are given in Table [T] 

The distribution on the sky of the stars included in the 
colour-magnitude diagram is illustrated in Fig. [4] 



4 MEMBERSHIP AND METALLICITY 

Babusiaux et al. (2005) have drawn attention to the con- 
siderable breadth in the infrared of Sculptor's giant branch. 
They deduce a range of 0.75 dex in [Fe/H] on the basis of 
fits of Padova theoretical isochrones for a fixed age to the 
Ks, V — Ks colour-magnitude diagram. Since then, a num- 
ber of kinematic and membership studies have been con- 
ducted, some of which also determine metallicities for indi- 
vidual stars. We use these new results to investigate the star 
formation history of Sculptor further. 

A proper-motion study covering 1177 stars within a 
7.25 arcmin radius of the Sculptor centre was conducted 
by Schweitzer et al. (1995). Our limiting magnitude is sig- 
nificantly brighter than theirs, but we have observed 407 
star in common with their presumably unbiased selection, 
including 381 member stars. 

Coleman et al. (2005) conducted a wide-area (10 square 
degrees centred on the galaxy) photometric and spectro- 
scopic survey of Sculptor with a view to mapping the outer 



structure of the galaxy. The V, V — I colour-magnitude dia- 
gram was used to select likely members. Velocities were de- 
termined from the spectra and the equivalent width of two 
Ca II triplet lines used to determine [Fe/H] on the Zinn & 
West (1984) metallicity scale. Data were only published for 
members, of which 84 are in common with our photometry. 

Another wide-area investigation was conducted by 
Westfall et al. (2006), this time employing Washington 
M, T2+DD051 photometry to isolate likely members. Their 
final catalogue contains 179 stars (157 velocity members) of 
which 113 are in common with our catalogue. No metallicity 
information was provided. 

Sculptor was included in a velocity survey of four Local 
Group dwarf spheroidals by Walker et al. (2009). The area 
covered extends well beyond our survey fields, and a total 
of 1541 stars with velocities and line strength information 
are catalogued. These were selected on the basis of their 
positions within a box drawn on a V, V — I colour- magnitude 
diagram. The bright limit appears to lie well below the tip 
of the red giant branch. Our survey has 762 stars in common 
with Walker et al., including 729 members. 

In a large-scale survey of abundances in Sculptor, Kirby 
et al. (2009) catalogue [Fe/H] values based on Ca II triplet 
measurements for 388 stars which were selected from the 
Westfall et al. (2006) paper. Only radial-velocity members 
are listed, but the underlying sample is subject to the same 
selection criteria as in the Westfall et al. study. Kirby et 
al. show that their abundances are in good agreement with 
those derived from high-resolution spectra by Battaglia et 
al. (2008). 

Finally, in a paper aimed at comparing Ca II triplet and 
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Figure 5. Colour-magnitude diagram for Sculptor with symbols 
as in Fig.[T] The line shows the expected position of the tip of the 
red giant branch assuming a distance modulus of 19.46 mag (see 
text). Definite proper motion or radial velocity members are over- 
plotted as green dots, while non-members are shown as magenta 
open squares. 



high-resolution spectroscopy of Sculptor stars, Battaglia et 
al. (2008) catalogue 93 stars, all of which are in common 
with our survey. They show that Ca II triplet abundances 
calibrated with respect to globular cluster [Fe/H] data give a 
reliable measure of metal abundance for members of Sculp- 
tor. 

Taking all of these surveys together, we find 952 definite 
members and 48 non-members in our sample, with a further 
796 stars of unknown status. Our final JHKs data, with 
membership status, are listed in Table The coordinates 
are as determined by us, and the column labelled "Scl#" 
refers to our internal numbers for the stars measured. 

The Ks , J — Ks colour- magnitude diagram is shown 
in Fig. O with the definite members marked in green, and 
definite non-members in magenta. Given the biases in the 
catalogues referred to above, too much should not be read 
into the distribution of stars with Ks magnitude, especially 
at the bright end where AGB stars, if any, have probably 
been discriminated against. The sloping line indicates where 
the tip of the red giant branch is expected to be (Bellazzini 
2008) for a range of J — Ks magnitudes as proxies for metal 
abundance. 

It is noteworthy that some of the 'definite' members 
lie some distance from the giant branch in Fig. [5] on the 
blue side. There are three outliers more than 0.1 mag bluer 
than the blue edge of the giant branch. Three stars from the 
Schweitzer et al (1995) catalogue, our numbers 1017, 2034 
and 9033, have probabilities of membership of 0.61, 0.90 and 
0.78, respectively. They all lie in a part of the J — H,H — K 
diagram that suggests they are most likely field stars, but 
it would be worth checking on this. On the red side, star 
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Figure 6. (a) Comparison of [Fe/H] for Battaglia et al. 
(2008) (dots) and Coleman et al. (2005) (triangles) with those of 
Kirby et al. (2009) for stars in common with our survey, (b) Sim- 
ilar comparison for Walker et al. (2009) ([Fe/H] from Mg indices) 
and Kirby et al. (2009) data. 



7082 is about 0.2 mag redder than the red edge of the giant 
branch, but is a member according to Walker et al. (2009). 

We have already alluded to the good agreement between 
the Kirby et al. (2009) and Battaglia et al. (2008) metallic- 
ities. The same applies to the Coleman et al. (2005) results, 
once they have been converted to the Carretta & Gratton 
(1997) scale. This is illustrated in Fig. |6]». Walker et al. 
(2009) measured equivalent widths for several Fe and Mg 
absorption lines in their spectra. They calibrated these with 
respect to globular cluster abundances to give [Fe/H] for the 
Sculptor stars. They cautioned against the use of these as 
absolute abundances since in two other Local Group galax- 
ies, abundances derived from Fe and Mg lines in this way 
gave rise to significantly narrower ranges than did those de- 
rived from the Ca II triplet. We find the same result for 
stars we have in common with both Kirby et al. (2009) and 
Walker et al. (2009) as shown in Fig.^b). 

We have determined the mean J — Ks as a function of 
Ks for the giant branch members and have divided them 
into two groups, redder or bluer than the mean. Using the 
metallicities from Coleman et al. (2005), Kirby et al. (2009) 
and Battaglia et al. (2008) we find a mean [Fe/H] of -1.48 
for 193 giants on the red side, and [Fe/H] of -1.81 for 184 
giants on the blue side. The two distributions are relatively 
broad and overlap, but the large range of abundances across 
the giant branch is clear. 

In Fig. [7] we show the colour-magnitude diagram for the 
members and C stars (MAG29 is off the red edge). Super- 
posed on the plot are isochrones (Marigo et al. 2008) for 
populations of age 12.7 Gyr and metal abundances, [Fe/H] 
= -1.76 (blue) and -1.4 (red), to illustrate that the width of 
the giant branch can be explained this way. Each isochrone 
has two parts, the rightmost one representing stars on the 
RGB, and the other the stars on the AGB. 

The presence of two very red Miras, both carbon stars, 
indicates that there must be an intermediate age population 
in Sculptor as well as the old population (see section 6). It is 
puzzling that the other carbon stars lie below the expected 
TRGB in the colour-magnitude diagram. However, Revaz et 
al. (2009) have found evidence for a population of age ~2 
Gyr. In Fig. [8] we show an isochrone for a 2 Gyr population 
and [Fe/H] ~ — 1.1 superposed on the colour-magnitude di- 
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Table 3. Distance modulus determinations for Sculptor. 

T 1 ' 1 ' I 




1.5 



J-K s 

Figure 7. Colour-magnitude diagram for Sculptor members. 
Symbols as in Fig. \E\ Isochrones for 12.7 Gyr populations with 
[Fe/H]=— 1.76 (blue lines) and —1.4 (red lines) are overplotted. 
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Figure 8. Colour-magnitude diagram for Sculptor members. 
Symbols as in Fig. \5\ An isochrone (Marigo et al. 2008) for a 
population of age 2 Gyr and [Fe/H] ~ —1.1 is overplotted as red 
lines. The region on the isochrone where C/O > 1 is shown in 
black. 



agram of members. The section of the isochrone marked in 
black shows where the C/O ratio exceed 1.0, which suggests 
that carbon stars can be produced in the region of the di- 
agram where they are observed if such a young relatively 
metal-rich population exists in Sculptor. 



Method 


Sculptor 


LMC 


554 d C Mira (M bol ) 


19.57 


(18.39) 


189 d C Mira (M bol /M K ) 


19.19 


(18.39) 


RR Lyracs (K) 


19.37 


18.37 


RR Lyraes (V) 


19.46 


18.38 


TRGB(V/) 


19.46 


(18.39) 



5 MIRAS AND THE DISTANCE TO 
SCULPTOR 

Table lists various estimates of the distance to Sculptor. 
The distances of the two carbon-rich Miras are based on 
the K and/or Mbol period-luminosity (PL) relations: 

M K = -3.52[logP - 2.38] - 7.24 (1) 
from Whitelock, Feast and van Leeuwen (2008), and 

M bol = -3.31\logP - 2.5] - 4.271 (2) 

from Whitelock et al. (2009). 

Both these relations were derived from carbon-rich Mi- 
ras in the LMC and assume a modulus of 18.39 mag obtained 
from a classical Cepheid distance scale based on parallaxes 
of Galactic Cepheids and with a metallicity correction ap- 
plied to the LMC Cepheid results (van Leeuwen et al. 2007). 
In the case of the 554-day Mira only the M BO i result is listed. 
In view of the redness of the star (J — K = 2.89) mag we an- 
ticipate significant circumstellar absorption at K and indeed 
the apparent distance modulus given using Eq. (2) is 0.4 mag 
greater than the Mbol modulus. In the case of the 189-day 
variable, Eq. (1) gives a modulus of 19.20 mag and Eq. (2) 
a modulus of 19.18 mag. Data for Sculptor RR Lyraes in K 
were taken from Pietrzyhski et al. (2008). The V data for the 
Sculptor RR Lyraes are from Kaluzny et al. (1995). A red- 
dening of E(B -V) = 0.018 mag from Schlegel et al. (1998) 
is adopted (see Pietrzyhski et al. 2008). The derivation of 
the PL(i^) and M(V)-[Fe/H] relations used is discussed by 
Feast (2010). The table also gives the LMC moduli obtained 
using the same relations. A recent discussion of the TRGB 
distance scale (Feast 2011) gives a distance modulus of 19.43 
mag for the LMC and 19.50 mag for Sculptor (using data 
from Sakai et al. (2000) and Rizzi et al. (2007)). Since we 
are here concerned primarily with the relative results of dif- 
ferent indicators we scale the Sculptor modulus in the table 
to an LMC value of 18.39 mag. 

The mean of the last three values is 19.43 mag and 
the mean of the two Miras is 19.38 mag, so the agreement is 
satisfactory. There is some concern that the 189-day variable 
gives a modulus 0.24 mag smaller than the mean of the non- 
Mira distance indicators. However, the scatter in the LMC 
PL relations is 0.15 mag for K and 0.12 mag in M 00 i (see 
Whitelock et al. 2008, 2009), so the modulus difference is 1.6 
to 2.0 a and may not be significant. One might suspect a 
difference could result from a metallicity effect on the zero- 
points, especially as the 189-day variable is likely to be very 
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Table 4. Metal abundances in Local Group galaxies. 



Galaxy [Fe/H] Mira Periods (days) 

Phoenix -1.37/-1.8 (old pop) 425 
Leo I -1.4 158, 180, 191, 252, 

283, 336, 523 

Fornax -0.9 (-2.0 to -0.4) 215, 258, 267, 280, 

350, 400, 470 

Sculptor -1.6 (14 Gyr) to -1.2 (4 Gyr) 189, 554 



old (see below) and metal poor. The LMC C-Mira relations 
are poorly denned at the short period end. However, short 
period Miras in Leo I and Fornax fit them well (Menzies et 
al. 2010, Whitelock et al. 2009). 



6 AGES OF CARBON-RICH MIRAS 

Mira variables have considerable potential both as Galactic 
and extragalactic distance indicators and as indicators of the 
age of stellar systems. Our work on dwarf spheroidals of the 
Local Group allows us to test Mira PL relations for possible 
metallicity dependence and to check whether our results are 
consistent with an period-age relation. 

For oxygen-rich Miras the evidence (summarized by 
Feast 2009) indicates that periods increase with decreasing 
age and increasing initial mass, from ~ 150 days at an age 
of ~ 12 Gyr to 450 days at ~ 3 Gyr and OH/IR Miras with 
periods ~ 1000 days are even younger. A broadly similar 
trend of period with age is found amongst the longer period 
carbon-rich Miras. Thus Galactic kinematics (Feast et al. 
2006) indicates an age of ~ 1.8 Gyr for carbon-rich Miras 
of period ~ 520 days with some indication of an increasing 
period with decreasing age. There is a Mira with a period 
of 450 days in NGC 1783 for which Milone et al. (2009) give 
an age of ~ 1.5 Gyr, and if a carbon-rich Mira of period 680 
days is a member of an LMC cluster (van Loon et al. 2003) 
its age is ~ 1.0 Gyr. It is likely that there will be a range 
in periods at a given age. How large this is is not precisely 
known. The cluster NGC 1978 (age ~ 2 Gyr (Milone et al. 

2009) ) contains two carbon-rich Miras for which Kamath et 
al. (2010) give periods of 376 and 458 days, while Nishida 
et al. give 491 days for the second star. In NGC 419 there 
are also two Miras (periods 488 and 738 days (Kamath et al 

2010) for the first of which Nishida et al. give a period of 526 
days) . The mean age of this cluster is ~ 1.5 Gyr though there 
appears to have been an extended period of star formation in 
it (Rubele, Kerber & Girardi 2010), which seems a common 
phenomenon in intermediate age clusters in the Magellanic 
Clouds (Milone et al. 2009). A more detailed discussion of 
the evolution of AGB variables will be given elsewhere, but 
it will be clear that Mira variables with periods ~ 200 days 
and lying on the PL relations with Mbol ~ — 3.6 mag are far 
below the RGB tip of clusters such as NGC1978 and must 
presumably have evolved from lower mass (older) stars. 

Table [4] lists the periods of the Miras in the four Local 
Group galaxies we have so far studied. All these Miras are 
known to be, or are very likely to be carbon rich. We can 
now ask how consistent are these Local Group results with 
the period-age correlation just described. 



Consider first the long period Miras (~ 400 — 550 days). 
These imply a significant population of a few Gyr age in 
all the systems. This is consistent with other evidence for 
Phoenix, Leo I and Fornax. 

In Phoenix there is evidence of significant star forma- 
tion in the age range 1 to 5 Gyr (see, e.g. Young et al. 2007, 
fig. 6) and a lower rate between 5 and 10 Gyr. In Leo I the 
star formation rate is relatively high in this same age range 
(e.g. Dolphin 2002, fig. 15). In Fornax, too, there is signif- 
icant star formation in this age range (e.g. Coleman & de 
Jong 2000, fig. 7). 

The case of Sculptor is particularly interesting in that it 
is frequently said to have only a very old population. Never- 
theless, the presence of a 554-day carbon-rich Mira indicates 
a significant population of age ~ 1 — 2 Gyr. Interestingly, 
the model of Revaz et al. (2009, fig. 13) does in fact show 
a small population in the 0-2 Gyr age range. Earlier work 
(Tolstoy et al. 2001, fig. 14) had this population at ~ 4 Gyr. 

In both Leo I and Fornax there has been significant star 
formation over the whole period, 1-12 Gyr. Thus, the range 
of Mira periods found in these galaxies is expected if period 
depends on age. This does not itself establish an age-period 
relation. However, as was previously pointed out (Menzies 
et al. 2010) two of the Leo I Miras fall in the outer parts 
of the galaxy where, as shown by Mateo et al. (2008), stars 
on the extended AGB are rare and where the main stellar 
population is older than in the inner parts. These two stars 
have periods of 191 and 283 days and are among the shorter 
period stars as we would expect for an older population. 

Sculptor is particularly interesting in this regard. Al- 
though, as we have just seen, the presence of a 554-day Mira 
is consistent with the small ~ 1 Gyr population of the model, 
there is no evidence of an intermediate population between 
that and the main population at ~ 10 Gyr and older. Since 
the evidence points against a carbon-rich Mira of 189 day 
period belonging to a young population we must assign it to 
this old population. 

Carbon stars in the Galactic halo are considered to be 
either intermediate age stars belonging to infalling streams 
or else old, low mass objects in binaries that have acquired 
their carbon from an initially more massive companion. This 
latter mechanism is reasonable since only a relatively low 
percentage of halo giants are carbon rich. This does not seem 
to be a viable explanation for the short period carbon-rich 
Miras in the dwarf spheroidals unless there are undetected 
short period oxygen-rich Miras in them, which is very un- 
likely. Evolutionary models of low mass stars do not usually 
produce carbon stars (e.g. Iben & Renzini 1983). The exact 
details of the mass and metallicity at which carbon stars do 
form depends on the models and we note that recent work 
(e.g. Karakas 2010; Suda & Fujimoto 2010) suggests that low 
mass stars with low metallicity may become carbon rich. 

It is also of interest to note that whereas the interme- 
diate age 554-day Mira is centrally located in Sculptor, the 
old 189 day one is outside the main concentration of stars 
in the galaxy. A similar Mira in Phoenix, with period 425 
days, is only about 1.3 arcmin from the centre of that galaxy, 
compared with the overall extent of about 8.7 arcmin (van 
de Rydt et al. 1991). 

The relative numbers of Miras in the various galaxies is 
roughly what we might expect from the total luminosities 
of the galaxies and estimates of the age distribution of their 
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stars. Thus, data in Amorisco & Evans (2010) suggest the 
total luminosities of Fornax and Leo I are 14±4 and 3.4± 1.1 
in units of 10 6 solar luminosities and both have had ongoing 
star formation from early times until quite recently. How- 
ever, the intermediate age population is relatively stronger 
compared to the very old population in Leo I than in For- 
nax (e.g. Dolphin 2002, fig. 15; Revaz 2009, fig. 13). Thus 
the roughly equal number of Miras is understandable. Sculp- 
tor is less luminous (~ 1.4 ± 0.6 in the same units,) while 
Phoenix has a similar luminosity, and one expects fewer Mi- 
ras in these galaxies. Both have a dominant old population 
and another increase in star formation later (< 5 Gyr in the 
case of Phoenix, < 2 Gyr in the case of Sculptor (Young et 
al. 2007; Revaz et al. 2009), though with the old population 
being relatively stronger in the case of Sculptor. Within the 
constraints of small number statistics these results for Mi- 
ras in Local Group galaxies are consistent with a general 
period-age relation for Miras and the known properties of 
the galaxies. 



7 CARBON-RICH MIRAS: METALLICITIES, 
MASS LOSS AND LUMINOSITIES 

The extent to which AGB stars produce dust at low metal- 
licities is of importance among other things for an under- 
standing of dust formation in the early universe. 

There is a general expectation that, other things being 
equal, the dust production in oxygen-rich AGB stars will de- 
crease with decreasing metallicity since the amount of avail- 
able material will decrease. In AGB carbon stars, however, 
this restriction does not necessarily apply since fresh carbon 
is transported to the stellar atmosphere by dredge-up. Tests 
of these expectations using observations in the Galaxy and 
the Magellanic Clouds have not been entirely conclusive (see 
Sloan et al. 2008; van Loon et al. 2008). Nevertheless, this 
and our other papers on dwarf spheroidals have shown the 
presence of carbon Miras with thick dust shells in metal-poor 
systems. Furthermore, Sloan et al. (2009) have used Spitzer 
Space Telescope spectra to deduce that Sculptor MAG29 is 
producing significant amounts of dust and have suggested 
on this basis that evolved low metallicity stars may make 
a significant contribution to dust production in the early 
Universe. 

Our present discovery that this star is a long period 
Mira allows us to further discuss this matter together with 
our work on Miras in other dwarf spheroidals. On the basis 
of the discussion of the last section, MAG29, with a period 
of 554 days, will be a member of the younger population in 
Sculptor. This has a mean metallicity of [Fe/H] ~ —1.4 (Tol- 
stoy et al.). This star is similar to the longest period Mira in 
Leo I, with a period of 523 days, and is thus of similar age. 
Gullieuszik et al. (2009) find from Ca II triplet observations 
of red giants that Leo I has a rather small dispersion in metal 
abundance, with [M/H] = -1.2, which they take as equiva- 
lent to [Fe/H] = -1.4 and thus similar to that inferred for 
MAG29. We have already shown (Menzies at al. 2010) that 
this star has a significant dust shell since, whilst it fits an 
Mboi-\og P relation, it is too faint, due to shell absorption, 
for a K-\og P relation. The difference between the apparent 
moduli for the K and M bo i relations gives an indication of 
the optical depth of the shell. For the Leo I star it is 0.8 mag 



whereas for MAG29 it is 0.4 mag indicating a thicker shell 
in the case of the Leo I Mira. There are other Miras in both 
Leo I and Fornax with similar thick shells indicating high 
mass-loss. Evidently MAG29 is not particularly exceptional 
in this regard. 

It is worth noting that in comparing mass-loss rates in 
different stellar systems it is important to consider similar 
types of objects. At least in the case of oxygen-rich Miras 
the mean mass-loss rate is a function of both luminosity (i.e. 
period) and pulsation amplitude (Whitelock et al. 1987). 

Groenewegen, Lancon & Marescaux (2009) have drawn 
attention to the great strength of the 1.53 /im band of C2H2 
(acetylene) band in MAG29. Sloan et al. have discussed the 
bands of this molecule at 7.5 and 13.7 /im in MAG29 and 
connect their strength with the low metallicity. The 1.53 
ixm band is known from the work of Joyce (1998) to vary 
strongly with pulsation phase in the 421-day carbon-rich 
Mira V Cyg and comparing his results with the AAVSO 
light curve shows it to be very strong near minimum light. 
Our infrared light curves show that MAG29 was near max- 
imum light when the data of both Groenewegen et al. and 
Sloan et al. were obtained implying unusual strength com- 
pared with stars so far observed in this region. However, 
we might expect it to be of similar strength in long-period, 
dust-enshrouded Miras in other Local Group systems. It is 
thus clear that carbon-rich Miras can produce significant 
amounts of dust in metal-poor systems. However, it is im- 
portant to realize that the carbon-rich Miras in the dwarf 
spheroidals we have discussed are, judging from their periods 
(< 500 days), too old (~ 1 Gyr) and of too low initial mass 
(~ 2Mq) to be directly relevant to the early Universe. At 
the metallicities of our Galaxy or the LMC, Miras with in- 
termediate initial masses and ages of relevance for the early 
Universe (periods > 1000 days) are not carbon rich. This 
is attributed to Hot Bottom Burning which turns the car- 
bon into nitrogen. Nevertheless, theory (e.g. Karakas 2010) 
suggests that at much lower metallicities, intermediate mass 
AGB stars do have C/O > 1 and will thus be carbon stars. 



8 CONCLUSIONS 

Infrared monitoring of Sculptor over a four year period has 
led to the following results: 

(1) A comparison with spectroscopically measured abun- 
dances shows that the width of the giant branch is primarily 
due to a metallicity spread with extreme metal poor stars 
([Fc/H] < -3.4) on the extreme blue side of the branch. 

(2) Apart from the two variables discovered, the known car- 
bon stars in the galaxy all lie below the TRGB. This con- 
trasts with the situation in other dwarf spheroidal galaxies. 
However, an isochrone with [Fe/H] w —1.1 and an age of 2 
Gyr (Marigo et al. 2008) has an abundance ratio C/O > 1 
extending down to the luminosities of these stars. 

(3) Two carbon-rich Miras are present in Sculptor and give 
a mean distance modulus of 19.38 mag based on an LMC 
modulus of 18.39 mag. 

(4) A general discussion of the relation of Mira periods to 
age, together with other evidence leads to the conclusion 
that the 554-day period Mira in Sculptor belongs to the 
small population with an age of ~ 2 Gyr whilst the 189-day 
Mira belongs to the dominant population of age ~ 10 Gyr. 
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(5) The 554-day Mira is identical with the star MAG29 
which has been previously discussed from Spitzer Space 
Telescope spectra as a low metallicity star with heavy mass 
loss. Using our previous work we find that there are sev- 
eral carbon-rich Miras in other metal-poor dwarf spheroidals 
with similarly thick circumstellar shells. Whilst this appears 
to strengthen the suggestion of Sloan et al. (2009) that low- 
metallicity carbon stars could make a significant contribu- 
tion to dust formation in the early Universe, this depends on 
whether young, high-mass AGB stars can become carbon- 
rich. 



9 ACKNOWLEDGEMENTS 

This research has made use of Aladin. This publication 
makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Mas- 
sachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the Na- 
tional Aeronautics and Space Administration and the Na- 
tional Science Foundation. This material is based upon 
work supported financially by the South African National 
Research Foundation. We are grateful to our colleagues, 
Toshihiko Tanabe, Yoshifusa Ita, Shogo Nishiyama, Enrico 
Olivier, Ryo Kandori and Jun Hashimoto, who assisted us 
or obtained some images of Sculptor for us; to Gary Da 
Costa and Giuseppina Battaglia for providing some useful 
data; and also to Greg Sloan for sending us the date of the 
Spitzer observation of MAG29. 

REFERENCES 

Amorisco N. C, Evans N. W., 2010, MNRAS, in press 
Azzopardi M., Lequeux J., Westerlund B. E., 1985, A&A, 144, 
388 

Azzopardi M., Lequeux J., Westerlund B. E., 1986, A&A, 161, 
232 

Baade W., Hubble E., 1939, PASP, 51, 40 

Babusiaux C, Gilmore G., Irwin M., 2005, MNRAS, 359, 985 

Battaglia G., Helmi A., Tolstoy E., Irwin M., Hill V., Jablonka 

P., 2008, ApJ, 681, 13 
Bellazzini M., 2008, Mem. S. A. It., 79, 440 
Carretta E., Gratton R. G., 1997, A&AS, 121, 95 
Coleman M. G., Da Costa G. S., Bland-Hawthorn J., 2005, AJ, 

130, 1066 

Coleman M. G., de Jong J. T. A., 2000, ApJ, 685, 933 

Cutri R. M., et al., 2003, University of Massachusetts and Infrared 

Processing and Analysis Center (IPAC/California Institute of 

Technology), 2MASS 
Dolphin A. E., 2002, MNRAS, 332, 91 
Feast M.W., 2000, J. Hist. Ast., 31, 29 

Feast M. W., Whitelock P. A., Menzies J. W., 2006, MNRAS, 
685, 933 

Feast M. W., 2009, in: (eds.) Toshiya Ueta, Noriyuki Matsunaga, 
Yoshifusa Ita, "AGB Stars and Related Phenomena, A Con- 
ference in Honor of Yoshikazu Nakada" , National Astronom- 
ical Observatory of Japan, Mitaka, Tokyo, Japan, p. 48 

Feast M. W., 2010, in: (eds.) Sterken C, Samus N., Szabados L., 
"Variable Stars, the Galactic halo and Galaxy Formation", 
Sternberg Astronomical Institute of Moscow University, Rus- 
sia, p. 45, larXiv:0912.4159l 

Feast M. W., 2011, in: ed. Gilmore G., "Planets, Stars and Stellar 
Systems", vol 5, Springer, Dordrecht, in press 



Frebel A., Kirby E. N., Simon J. D., 2010, Nature, 464, 72 
Frogel J. A., Blanco V. M., McCarthy M. F., Cohen J. G., 1982, 
AJ, 252, 133 

Groenewegen M. A. T., Lancon A., Marescaux M., 2009, A&A, 
504, 1031 

Gullieuszik M., Held E. V., Saviane I., Rizzi L., 2009, A&A, 500, 
735 

Hurley-Keller D. A., 2000, PhD thesis, University of Michigan 
Iben I., Renzini A., 1983, ARAA, 21, 271 
Joyce R. R., 1998, AJ, 115, 2059 

Kaluzny J., Kubiak M., Szymahski M., Udalski A., Krzcminski 
W., Mateo M., 1995, A&AS, 112, 407 

Kamath D., Wood P. R., Soszyfiski I., Lebzelter T., 2010, MN- 
RAS, 408, 522 

Karakas A. I., 2010, MNRAS, 403, 1413 

Kato D., ct al. 2007, PASJ, 59, 615 

Kirby E. N., Guhathakurta P., Bolte M., Sneden C, Geha M., 

2009, ApJ, 705, 328 

Marigo P., Girardi L., Bressan A., Groenewegen M. A. T., Silva 

L., Granato G. L., 2008, A&A, 482, 883 
Mauron N., Azzopardi M., Gigoyan K., Kendall T. R., 2004, 

A&A, 418, 77 

Mateo M., Olszewski E. W., Walker M. G., 2008, ApJ, 675, 201 
Menzies J., Feast M., Tanabe T., Whitelock P., Nakada Y., 2002, 

MNRAS, 335, 923 
Menzies J., Feast M., Whitelock P., Olivier E., Matsunaga N., da 

Costa G., 2008, MNRAS, 385, 1045 
Menzies J. W., Whitelock P. A., Feast M. W., Matsunaga N., 

2010, MNRAS, 406, 86 

Milone A. P., Bedin L. R., Piotto G., Anderson J., 2009, A&A, 
497, 755 

Pietrzyriski G., Gieren W., Szewczyk O., Walker A., Rizzi L., 
Bresolin F., Kudritzki R-P., Nalewajko K., Storm J., Dall'Ora 
M., Ivanov V., 2008, AJ, 135, 1993 

Revaz Y., Jablonka P., Sawala T., Hill V., Letarte B., Irwin M., 
Battaglia G., Helmi A., Shetrone M. D., Tolstoy E., Venn K. 

A. , 2009, A&A, 501, 189 

Richer H. B., Westerlund B. E., 1983, ApJ, 264, 114 

Rizzi L., Tully R. B., Makarov D., Makarova L., Dolphin A. E., 

Sakai S., Shaya E. J., 2007, ApJ, 661, 815 
Rubele S., Kerber L., Girardi L., 2010, MNRAS, 403, 1156 
Sakai S., Zaritsky D., Kennicutt R. C. Jr., 2000, AJ, 119, 1197 
Salaris M., Girardi L., 2005, MNRAS, 357, 669 
Schechter P. L., Mateo M., Saha A., 1993, PASP, 105, 1342 
Schlegel A. J., Finkbeincr D. P., Davies M., 1998, ApJ, 500, 525 
Schweitzer A. E., Cudworth K. M., Majewski S. R., Suntzeff N. 

B. , 1995, AJ, 110, 2747 
Shapley H., 1938, Harv. Bull., 908 

Sloan G. C, Kracmer K. E., Wood P. R., Zijlstra A. A., Bernard- 
Salas J., Devost D., Houck J. R., 2008, ApJ, 686, 1056 

Sloan G. C, Matsuura M., Zijlstra A. A., Lagadec E., Groenewe- 
gen M. A. T., Wood P. R., Szyszka C, Bernard-Salas J., van 
Loon J. Th., 2009, Science, 323, 353 

Suda T., Fujimoto M. Y., 2010. larXiv:1002.0863l 

Thackeray A.D., 1950, Obs., 70, 144 

Tolstoy E., Irwin M. J., Cole A. A., Pasquini L., Gilmozzi R., 

Gallagher J. S., 2001, MNRAS, 327, 918 
Tolstoy E., Irwin M. J., Helmi A., Battaglia G., et al., 2004, ApJ, 

617, L119 

van Agt S. L. Th. J., 1978, Publ. David Dunlap Obs., 3, 205 
van der Rydt F., Demers S., Kunkel W. E., 1991, AJ, 103, 130 
van Leeuwen F., Feast M. W., Whitelock P. A., Laney C. D., 

2007, MNRAS, 379, 723 
van Loon J. Th., Marshall J. R., Matsuura M., Zijlstra A. A., 

2003, MNRAS, 341, 1205 
van Loon J. Th., Cohen M., Oliveira J. M., Matsuura M., McDon- 
ald I., Sloan G. C, Wood P. R., Zijlstra A. A., 2008, A&A, 
487, 1055 



AGB Stars in the Sculptor Dwarf Spheroidal Galaxy 



Walker M. G., Mateo M., Olszewski E. W., 2009, AJ, 137, 
Westfall K. B., Majewski S. R., Ostheimer J. C, Frinchaboy P. 

M., Kundel W. E., Patterson R. J., Link R., 2006, AJ, 131, 3 
Whitelock P. A., Feast M. W., Pottasch S. R., 1987, in: (eds.), 

"Late stages of stellar evolution" , Reidel, Dordrecht, p. 269 
Whitelock, P. A., Feast, M. W., Marang, F., Groenewegen, M. A. 

T., 2006, MNRAS, 369, 751 
Whitelock P. A., Feast M. W., van Leeuwen F., 2008, MNRAS, 

386, 313 

Whitelock P. A., Menzies J. W., Feast M. W., Matsunaga N., 

Tanabe T., Ita Y., 2009, MNRAS, 394, 795 
Young L. M., Skillman E. D., Weisz D. R., Dolphin A. E., 2007, 

ApJ, 659, 331 
Zinn R., West M. I., 1984, ApJS, 55, 45 



